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Abstract. Based on a new and large sample of op- 
tical and near-infrared spectra obtained at the Mount 
Stromlo and Siding Spring Observatories (Langon & Wood 
1998; Langon & Wood, in preparation), spectrophotomet- 
ric properties of cool oxygen- and carbon-rich Long Pe- 
riod Variables and supergiants are presented. Tempera- 
tures of oxygen-rich stars are assigned by comparison with 
synthetic spectra computed from up-to-date oxygen-rich 
model atmosphere grids. The existence of reliable optical 
and near-infrared temperature indicators is investigated. 
A narrow relation between the bolometric correction BCi 
and the broad-band colour I— J is obtained for oxygen- 
rich cool stars. The ability of specific near-infrared indices 
to separate luminosity classes, atmospheric chemistry or 
variability subtypes is discussed. Some comments are also 
given on extinction effects, water band strengths in Long 
Period Variables and the evaluation of 12 CO/ 13 CO ratio 
in red giants. 

Key words: Stars: AGB and post-AGB - Stars: carbon - 
Stars: fundamental parameters - Stars: late-type - Stars: 
oscillations 



1. Introduction 

The near-infrared (NIR) radiation of our and other galax- 
ies is essentially provided by luminous red stars: red super- 
giants dominate in starburst regions, Asymptotic Giant 
Branch stars (AGBs) , among which the O-rich and C-rich 
Long Period Variables (LPVs) of various variability types, 
take over in stellar populations with ages between 10 8 and 
10 9 years, and First Giant Branch stars dominate in older 
regions (Bruzual & Chariot 1993; Langon 1999). In view 
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of the recent developments in NIR observational technol- 
ogy, tools have become necessary that would allow us to 
connect NIR spectrophotometric signatures with the fun- 
damental stellar properties provided by theoretical stellar 
evolution tracks. Only then can integrated NIR emission 
properties or star counts be interpreted in terms of galaxy 
history. These connection tools are also needed from a 
purely stellar perspective. The evolution, structure and at- 
mospheres of the evolved red stars are among the standing- 
problems in stellar astrophysics and valuable constraints 
on the models are to be expected from population stud- 
ies. The choice of suitable filter passbands, able to separate 
stellar subtypes, is a prerequisite for efficient surveys. 

In this paper, the emphasis is set on LPVs, which con- 
stitute the most luminous and coolest but also the most 
complex population of red giants in regions where more 
massive red supergiants have died. Static red giant model 
atmospheres are progressing rapidly due to recent im- 
provements in molecular opacity data (Plez 1999a, 1999b). 
Variable star models are still far from satisfactory because 
of the interplay between pulsation, radiative transfer and 
chemistry (see Hofner 1999 for a review). Empirical spec- 
troscopic or spectrophotometric studies adequate for our 
purpose have remained scarce because they are time con- 
suming and because a significant wavelength coverage re- 
quires synchronised observations: LPV periods are typi- 
cally one year long and cycle to cycle variations common. 
In addition, since good infrared instrumentation is usually 
found on relatively large telescopes that tend to be pref- 
erentially directed towards faint sources instead of bright 
stars, most previous photometric studies of LPVs have 
been optical (Kerschbaum & Hron 1996; Alvarez & Men- 
nessier 1997; Wing et al. 1998). Based on a new, rich sam- 
ple of complete optical+NIR spectra of red stars, we are 
able to extend these studies. 
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The properties of the data are summarised in Sect. 2. 
In Sect. 3, we discuss spectrophotometric temperature in- 
dicators at optical and NIR wavelengths. The physical sig- 
nificance and absolute calibration of the empirical scales, 
in particular for variable stars, will be discussed in more 
depth elsewhere. The subsequent sections discuss the abil- 
ity of commonly used and newly defined NIR narrow band 
indices to separate luminosity, surface chemistry or vari- 
ability subtypes. The results are summarised in Sect. 10. 

2. The sample 

The reader is refered to Lancpn & Wood (1998) and 
Langon & Wood, in preparation, for a thorough descrip- 
tion of the sample. Around one hundred spectra of cool 
stars ranging from 5100 A to 2.5 /jm have been obtained 
at the Mount Stromlo and Siding Spring Observatories, 
using the three cross-dispersed grisms of CASPIRj] on the 
2.3 m Telescope (near- infrared spectra with a resolution 
R of 1100) and a low resolution grating spectrograph on 
the 74 inch Telescope (optical spectra, R ~ 30) . The sam- 
ple mainly contains instantaneous spectra of O-rich and 
C-rich local LPVsf| but also includes local non- variable gi- 
ants and a few luminous LPVs from the Magellanic Clouds 
and the Galactic Bulge. Additional spectra restricted to 
either the optical (5100-9900 A) or the NIR region (9750- 
24900 A) are available for some stars and are included in 
the analysis whenever possible. 

3. Temperature indicators 

3.1. Temperature definitions 

Only little is known about effective temperatures of LPVs. 
In principle, the effective temperature relates the bolomet- 
ric luminosity of a star to its effective radius and is thus 
accessible through surface brightness measurements, i.e. 
apparent luminosities combined with angular diameters 
(e.g. Perrin et al. 1998). However, due to their very ex- 
tended atmospheres, LPVs such as Miras cannot easily be 
assigned a unique radius: observational determinations are 
wavelength dependent and depend on the adopted limb- 
darkening law (Labeyrie et al. 1977; Baschek et al. 1991; 
Hofmann et al. 1998). The use of model atmospheres is 
inevitable not only in providing these laws in the various 
observational filter passbands, but also to relate the re- 
sulting "monochromatic" radii to the effective radius (usu- 
ally based on Rosseland mean optical depths) that is used 
in the effective temperature definition of stellar models. 
Large uncertainties and inconsistencies between authors 
occur as a result of e.g. incomplete treatment of molecu- 
lar opacities in the models, of different assumptions on the 
stellar pulsation mode or on whether or not the phase in 

1 Cryogenic Array Spectrometer and Imager, McGregor 1994 

2 All variability types used in this paper are from the General 
Catalogue of Variable Stars (Kholopov et al. 1988) 



the pulsation cycle is taken into account (Scholz & Takeda 
1987; Haniff et al. 1995; van Belle et al. 1996; Bessell et 
al. 1996). 

Photometric indices, i.e. colours and molecular band 
strength measurements, provide alternative effective tem- 
perature estimates (Ramirez et al. 1997; Alvarez & Men- 
nessier 1997). Optical indices define the spectral types, 
which are calibrated against the effective temperature ei- 
ther on the basis of angular diameter measurements (Ridg- 
way et al. 1980; Perrin et al. 1998) or by direct comparison 
with synthetic spectra (Fluks et al. 1994). The calibrations 
inherit the model dependence already mentioned. In ad- 
dition they generally assume a unique relation between 
spectral type and effective temperature, independent of 
variability type or time. This assumption oversimplifies 
the real situation: Alvarez & Plez (1998) have shown that 
TiO and VO bands, two temperature sensitive indicators, 
do not vary in phase along the Mira luminosity cycle, and 
the data used here contains pairs of spectra with insignif- 
icant differences below 1 /im but very different NIR H2O 
vapour bands. As expected from the complex stratification 
of dynamical atmosphere models (Bowen 1988; Fleischer 
et al. 1992; Hofner et al. 1998) molecular indices contain 
local information and are sensitive each in its own way to 
the details of the pulsation history as well as to the global 
property called effective temperature. 

In view of the above difficulties, this paper avoids di- 
rect comparison with dynamical model atmospheres. A 
temperature refered to as the "static" effective temper- 
ature (and noted T e g throughout the paper for simplic- 
ity)]^] is nevertheless tentatively assigned to each spectrum, 
based on the fact that the optical part of the static and 
variable giant spectra in our sample can indeed be rea- 
sonably well fitted with static giant models. The resulting 
absolute T e ff scales may be biased estimates of the real 
effective temperatures, but a better understanding of pul- 
sating models is required before these biases can be eval- 
uated and their dependence on the pulsation properties 
accounted for. The static T e g scale is used to calibrate 
photometric temperature scales and all further discussion 
then relies on purely empirical findings. 

3.2. Static effective temperature assignment 

The synthetic spectra used for the static T ff assignment 
are computed from a large grid of MARCS models (Plez 
et al. 1992; Plez 1992; Plez, private communication) repre- 
senting red giants with solar composition and different val- 
ues of atmospheric extension, surface gravities log g rang- 
ing from —0.5 to 3.5, and of effective temperatures from 
2400 to 4750 K. The spectral synthesis includes the up- 
dated opacity data for TiO (Plez 1998), VO and H 2 as 
detailed in Bessell et al. (1998) and Alvarez & Plez (1998). 

3 "Effective temperature" , in extenso, is used when the the- 
oretical quantity is meant. 
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Fig. 1. Observed spectrum of SPhe, a M6 SRb vari- 
able (full line) and synthetic spectrum with T e ff=3400 K, 
log 0=1.0 and M.—1.0 Mq (dashed line). The dotted lines 
are two synthetic spectra with T e ff=3200 K and 3600 K 



The latest models are preliminary computations of a large 
new and updated grid now being computed by the Upp- 
sala group. 

The static T e g associated with a spectrum is defined 
here as the effective temperature of the MARCS model 
that provides the minimum \ 2 fit to the 5100-9950 A 
data, logg and atmospheric extension being allowed to 
vary freely. 

We have determined this T c g for 131 spectra of oxygen- 
rich giants, among which some are Mira, Semi-Regular 
(SRa and SRb types) or Irregular variables, observed at 
various phases. The sensitivity of the spectra to T c ff com- 
bined with the step of the model grid allows us to assign 
T c ff values to ±100 K (cf. Fig. [l]). The corresponding ac- 
ceptable logg values lie between 0.0 and 1.5, a range over 
which spectral variations are much smaller than those pro- 
duced by a 100 K change in T e g. 

3.3. Photometric temperature indicators 

The following sections discuss narrow band and broad 
band colours derived from the observed spectra. We use 
filters similar to those of Wing (1967) and White & Wing 
(1978), adopting the passbands defined by Bessell et al. 
(1989a); we complete the set with filters from Fluks et 
al. (1994) and additional indices designed to measure OH, 
H2O and CO band strengths. The VRIJHK colours are 
computed with the passbands of Bessell & Brett (1988) 
and Bessell (1990). The zero point calibrations are based 
on a spectrum of Vega (Castelli & Kurucz 1994) as used 
by Bessell et al. (1998): the colours of Vega are set equal 
to zero, except for the Wing colours for which we used 
Vega's observed colours (Wing, private communication). 



3.3.1. Optical temperature indicators 

Most optical colours appear to be quite good temperature 
indicators (e.g. 710-888, 847-839, 710-745, 888-883 as 
used by Wing, or Siy2,s p , S 2 /z,Sp as m Fluks et al. 1994). 

Figure ^ shows the 710—888 index as a function of the 
T e ff assigned previously. This index has been chosen to 
be our temperature indicator subsequently as its scale of 
variation is large and the scatter is small. It measures the 
ratio between the bandheads of the TiO 7 Av=0 and the 
TiO S Av=0 ro-vibrational transitions (Brett 1990). The 
710 and 888 filters are displayed on Fig. |[ A least-square 
polynomial fit gives: 



(710 - 888) = 5.8 + 3.2 10~ 4 x T cff 



4.7 10~ 7 x TL 



(1) 



The standard deviation is 0.3 mag in the range 2500 K < 
T cf f < 3900 K. 

It is worth noticing that non- variable and variable stars 
on one hand, local, Bulge and LMC/SMC stars on the 
other hand, all seem to follow the same relation. This 
may be surprising at first thought. In fact, it shows that 
the relation between the 710—888 index and the whole 
optical spectrum used in our \ 2 adjustment doesn't de- 
pend strongly on metallicity or variability type. The opti- 
cal parts of the spectra all tend to lie on a tight sequence 
(which makes the fits with static model spectra possible). 
It doesn't imply, however, that the effective temperature 
scale along this sequence is independent of the object type 
or metallicity: our T c g assignment based on the optical 
spectrum is liable to be a good estimate of the real effec- 
tive temperature of static stars, but could be systemati- 
cally different from the effective temperature in variable 
objects. 

Only two Bulge stars of our sample were observed 
in the optical region. More observations are required to 
make firm statements on their specific static T e g- versus 
(710-888) relation. 

The outlying points in Fig. | at 710-888 > 4 are 
three observations of SVLib, a M8 Mira star. These val- 
ues of the TiO index would lead to particularly low tem- 
peratures when extrapolating the T e g-(710— 888) relation 
without caution. However, the global energy distribution 
of the star and the discussion in the following section ar- 
gue against such an interpretation. These points were not 
taken into account when deriving Eq. ([!]). 

3.3.2. Near-infrared temperature indicators 

For extragalactic study purposes, near-infrared tempera- 
ture indicators are of great interest. Indeed, the contribu- 
tion of AGB stars to the emission of a stellar population is 
largest in this spectral range, and the population synthesis 
predictions are thus particularly sensitive to temperature- 
colour-spectral feature relations (Lancon 1999). 

We have applied various near-infrared filters to the ob- 
served spectra. In Fig. [|, the Wing VO index 105—104 is 
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Fig. 2. The 710—888 index as a function of the effective 
temperature derived from the fit of MARCS models to the 
optical spectrum (plus signs: M Miras; empty circles: M 
Semi-Regulars; small empty circles: M Irregulars; squares: 
K and M giants; triangles: supergiants; diamonds: bulge 
stars; crosses (X): LMC/SMC stars 
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Fig. 3. Narrow-band filters used in Sect. 3. The spectra 
are from bottom to top: RPhe, a O-rich Mira; RUPup, 
a C-rich Semi-Regular and HD 101712, a supergiant. Re- 
gions around 1.4 and 1.9 /im are associated with extreme 
telluric water absorption 



plotted against our optical T c s indicator 710—888. There 
is a good correlation between the two colours, the VO in- 
dex becoming a useful T c g index below 3300 K, as already 
found by Bessell et al. (1989a) from synthetic spectra. 
The 104 and 105 filters are shown in Fig. [| The scatter 
of Fig. ^| might be due to atmospheric extension, surface 
gravity or metal abundance (Bessell et al. 1989a) but the 
main source of scatter is certainly variability, Miras exe- 
cuting loops in such TiO /VO diagrams (Bessell et al. 1996; 
Alvarez & Plez 1998). 



Fig. 4. The 105—104 index versus the T c s indicator 
710-888 



The three outlying points (710— 888>4) are, again, the 
observations of SVLib. This M8 Mira star has strong TiO 
absorption at 710 nm but normal VO bands (see also 
Lancon & Wood 1998). The observations were taken at 
different dates: the effect cannot be associated with errors 
in the acquisition or the reduction of the data. This star 
deserves further scrutiny and in particular a more com- 
plete phase coverage to clarify its pecularities. 

It is difficult to define near-infrared colours more con- 
venient than 105—104 that might be used as T c s indi- 
cators. The CO (2.3 /itm) equivalent width suggested by 
Ramirez et al. (1997) appears inadequate: the correlation 
with T e ff is only marginally significant in our LPV sample. 
The large dispersion is mainly due to variations of the or- 
der of 10 A in the CO equivalent width of individual stars, 
which are not in phase with the T c g variations. Differences 
in metallicity and C/O ratios are likely to contribute to 
the dispersion. Figure || shows the indices most tightly 
correlated with 710-888. The 200-220 index which is a 
measure of the strength of H2O at 2.0 /im relative to the 
continuum point at 2.2 fim reveals a more pronounced 
scatter than the VO index when plotted against the T e s 
indicator 710-888. The 200 and 220 filters are displayed 
on Fig. § 

The broad-band colour V— K is not a pure near- 
infrared index. It is a very good effective temperature in- 
dex for static cool stars, not so much affected by metallic- 
ity effects according to Bessell et al. (1998), and we find 
it very well correlated with the 710—888 index. V— K val- 
ues have been computed after extrapolating each observed 
spectrum (which begins at 5100 A and ends at 24900 A) 
on the blue and on the red sides by the synthetic spectrum 
which gave the best x 2 fit. The correction is always less 
than 0.15 mag (mean: 0.09 mag). Bessell et al. (1998) have 
determined polynomial fits to empirical effective temper- 
ature versus V— K relations for giants. Figure || shows the 
comparison between the effective temperature obtained 
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Fig. 6. The effective temperature obtained from V— K and 
the relation of Bessell et al. (1998) against our static T e g 
assignment derived from the fit of MARCS models to the 
optical spectrum (see Sect. 3.2). The dotted line is the 
one-to-one line. Same symbols as Fig. [| 

using their relation (their Table 7, data sets 'abed') and 
the static T g we determined in Sect. 3.2: the agreement 
is good. The broad-band colours I— J and I— K are almost 
as good T e ff indicators as V— K (Fig. |^). A systematic dif- 
ference between the supergiants (triangles) and the other 
stars is seen: their I— J and 1— K tend to be larger at given 
710-888. 

4. Bolometric corrections 

The large spectral region covered by the observations and 
the fact that the cool stars emit most of their energy in 
the near-infrared enable us to compute bolometric cor- 
rections. A first order correction to the determination of 
the bolometric luminosity has been applied by extrapolat- 
ing the observed energy distributions on the blue and red 
sides by synthetic spectra as explained in Sect. 3.3.2. The 
mean correction on the bolometric magnitudes amounts 
to 0.22 mag. No additional correction was added in the 
case of dust shell enshrouded stars (see Sect. 5). The zero- 
point of the bolometric scale is computed with the follow- 
ing values (Bessell et al. 1998): M bo i,0=4.74, M v ,0=4.81, 
V Q =-26.76, f to t,0=l-371xlO 6 erg cm" 2 s" 1 . The abso- 
lute flux at 5556 A of the Vega model is set to 3.44 10 -9 
erg cm -2 s _1 A -1 (Hayes 1985) to achieve a V magnitude 
of 0.03 mag for Vega (Bessell et al. 1998). Figure @ presents 
the bolometric correction BCi = mboi~ I versus I— J. This 
relation is of particular interest for the ongoing infrared 
sky surveys like DENIS (Epchtein 1998). It is remarkable 
that local, Bulge and LMC/SMC stars, and all variabil- 
ity types follow the very same relation: there is a single 
narrow sequence for oxygen-rich stars. A polynomial fit 
gives: 

m b oi - I = 1-32 - 0.574 x (I - J) - 0.0646 x (I - J) 2 (2) 



Fig. 7. The bolometric correction BCi as a function of the 
I— J colour for cool oxygen-rich stars. The dotted line is 
the polynomial fit 

The standard deviation of the relation is 0.09 mag. 

A complete study of the bolometric corrections derived 
from this sample will be presented elsewhere (Mouhcinc 
et al., in preparation). 

5. Extinction 

The colour indices used in this paper have not been cor- 
rected for extinction. To what extent may interstellar and 
circumstellar matter bias the two-colour relations and add 
to their dispersion? 

In order to estimate the effects of interstellar extinction 
on the colours, we used the extended mean extinction law 
given by Fluks et al. (1994) based on the observed mean 
extinction law of Savage & Mathis (1979) and on the the- 
oretical extinction law of Steenman & The (1989, 1991). 
From the tabulated R(X) = a(X) / Eb-v values, we com- 
pute the interstellar extinction correction for band Bi as 
a function of Eb-v'- 



Bi — B 



i.O 



-2.5 log 



r °° s i (x)F(\)dx 



/„ S , l (A)F(A)10°- 4fi ( A ) £; «-vdA 



(3) 



where Bi is the uncorrected index, Bi$ is the dereddened 
index, Si(X) is the response function within the band Bi 
and F(X) is the observed spectrum. As an approximation, 
we take F(X) as the arithmetic mean of all the spectra. 

To first order, reddening mimics a shift towards lower 
temperatures in all the usual broad band two-colour plots 
(e.g. Fig. ||). But molecular band indices may be little 
affected or even vary orthogonally to the temperature se- 
quence, depending on the separation between the two pho- 
tometric passbands used, and on the side of the molecular 
absorption feature on which the "continuum" is measured. 
Thus, the effect of dust is expected to be most obvious in 
plots that combine a broad band colour and a molecular 
index. An example is given in Fig. |^. Most of the local 
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Fig. 5. The 200—220, V— K, I— J and I— K indices versus the temperature indicator 710- 
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Miras are grouped along a sequence at J— H ~ 0.85. The 
dispersion across the sequence is real, i.e. much larger than 
observational errors. But only a small fraction of the dis- 
persion may be associated with interstellar extinction: the 
motion of individual stars in J— H during their pulsation 
cycle is comparable with the width of the sequence (see the 
example of the M3 Mira X Men on Fig. ||), while the effect 
of typical values of the optical depth towards local stars 
is much smaller (comparable in amplitude to the observa- 
tional uncertainties). This justifies neglecting interstellar 
extinction towards local stars. 

Two local Miras with very large J— H values (larger 
than 1.5) are identified in Fig. ||. IO Vir, a M6 Mira, was 
observed twice. It is a known OH maser (OH 334.7+50.0, 
Turner 1979). This star is thus enshrouded in an im- 
portant dust shell. So is the second extreme Mira, 
WXPsc (OH 128.6-50.1). IOVir and WXPsc are the 
only two Miras of our sample with an OH designation: 
this corroborates the fact that circumstellar matter, as 
opposed to interstellar matter, is responsible for their 



red colours. Infrared excess is measured by the ratio 
F(12 /xm)/F(2.2 fixa) (see Habing 1996 and references 
therein; Le Bertre & Winters 1998). This ratio has been 
derived from the IRAS PSC (IRAS Science Team 1988) 
and from the K magnitude (Gezari et al. 1996), whenever 
available (using the calibration of Beckwith et al. (1976) to 
convert K magnitudes into fluxes at 2.2 /jm). As expected, 
Fig. [)] clearly shows that WX Psc and IO Vir (and to some 
extent IRC -20427) have large F(12 /im)/F(2.2 fim) ra- 
tio. The presence of cool circumstellar material around 
these stars is thus confirmed and explains their large J— H 
colours. 

Even in the extremely obscured cases mentioned 
above, the data do not indicate a significant contribu- 
tion of circumstellar continuum emission at wavelengths 
shorter than 2.4 /xm. Such a warm dust contribution would 
reduce the equivalent widths of the CO bands (Tanaka et 
al. 1996). No significant reduction is observed in our es- 
sentially O-rich sample. The depths and shapes of the CO 
bands (A > 2.29 /im) of the naturally reddened spectra 
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compare very well to those of cool unobscured Miras like 
XMen, once the latter have been artificially reddened to 
match the colours of the former. 

The Bulge stars systematically lie above the sequence 
of local optical LPVs. The interstellar extinction effect is 
obvious. The extinction towards the SgrI and NGC 6522 
galactic Bulge fields, to which the Bulge stars of our sam- 
ple belong, corresponds to Eb-v — 0.55 according to es- 
timates in the literature (Wood & Bessell 1983; Glass et 
al. 1995; Barbuy et al. 1998). Once corrected for redden- 
ing, the observed Bulge stars lie on the edge of the local 
Mira sequence. This is consistent with the fact that, for 
Bulge stars, a large J— K index has been used as a se- 
lection criterion for inclusion in the spectroscopic sample. 
However, the spectra indicate that at least some of the 
observed Bulge stars are reddened more than the aver- 
age SgrI and NGC 6522 field stars: for instance, the spec- 
trum of SgrI N°ll (observed in July 1996) is very well 
matched, in terms of the energy distribution, the molec- 
ular bands and the global aspect of metal line blends, by 
the spectrum of the local Mira RS Hya (February 1995) if 
Eb-v — 1- An equally good match is obtained between 
the local XMen (March 1996) and NGC 6522 N°435 (May 
1996) with Eb-v = f -6. This additional extinction may be 
due in various proportions to circumstellar matter and to 
the patchy distribution of interstellar dust. Similarly high 
extinction values are found towards various star clusters in 
the neighbourhood of the Galactic Centre (Barbuy et al. 
1998). Only two of the Bulge stars of the sample have been 
detected by IRAS (NGC 6522 N°435= IRAS 17593-3006, 
NGC 6522 N°205=IRAS 17591-2959; Glass 1986; Glass et 
al. 1995). Using the K magnitudes of Wood & Bessell 
(1983), we obtain 12 /im to 2.2 /im flux ratios of 1.05 and 
1.6, respectively, i.e. relatively low values as compared to 
the local IRAS sources of Fig. ||. Again, the strong CO 
absorption longward of 2.29 /zm argues against a signifi- 
cant contribution of warm circumstellar dust emission to 
the near-IR light of the Bulge LPVs. 

The Bulge stars tend to display strong VO and TiO 
bands; molecular bands need to be observed in a larger 
unbiased sample in order to determine whether this is or 
is not due to the selection of red and luminous objects. 
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Fig. 8. J— H versus the near-infrared T c g indicator 
105—104. The arrow gives the extinction vector for an 
Eb-v equal to 1. The cross in the upper left corner 
gives the typical error bars. Lines connect observations 
of X Men, an M3 Mira, at various phases. Same symbols 
as Fig. H 
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Fig. 9. The F(12 /zm)/F(2.2 fim) ratio versus J-H. Three 
stars (WXPsc, IO Vir and IRC -20427) present an in- 
frared excess due to the presence of a cool circumstellar 
shell. Same symbols as Fig. [| 



6. Water in Long Period Variable stars 

The most distinctive features in the near-IR spectra of 
pulsating red giants are their impressive H2O absorption 
bands. Johnson & Mendez (1970) already observed Mira 
type LPVs with much deeper water vapour bands than 
static giants of similar spectral types. LPV models (Bessell 
et al. 1989b; Bessell et al. 1996) and observations (e.g. 
Matsuura et al. 1999) have subsequently confirmed that 
deep H2O features are related to pulsation and have pro- 
vided basic understanding: pulsation produces much more 
extended atmospheres, and in addition dense cool layers 
may result from the periodical outwards running shocks. 



In various ways pulsation thus leads to the existence of 
regions where relatively low temperatures (~ 10 3 K or 
lower) are combined with relatively high densities, condi- 
tions that favour the formation of H2O and possibly even 
of dust (Bowcn 1988; Feuchtinger et al. 1993; Jeong et al. 
1999). 

Figure H (top left) shows that H2O absorption indices 
for LPVs range from to values far above those the coolest 
static giants reach. At a given stellar energy distribution a 
large range of H2O absorption indices is found (a similar 
phenomenon has been described for SiO bands by Aringer 
et al. 1999). Even relatively warm LPV spectra (instan- 
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taneous T c ft > 3500 K) may display deep water vapour 
bands, a behaviour that remains to be reproduced quan- 
titatively by LPV models (Mouhcine & Langon 1998). 
More data (in particular complete phase coverage for more 
stars) is required to gain better insight into the relation 
between the H2O features, pulsation amplitude, phase, 
metallicity and details of the stellar structure. 

For the purposes of this paper, we conclude that deep 
H2O features can be taken as a clear indication for pulsa- 
tion, while the absence of water vapour does not exclude 
that the observed star be an LPV. 

7. Supergiants, giants and variability types 

From a spectrophotometric point of view, supergiants, gi- 
ants and dwarfs are usually separated mainly on the basis 
of CO and H2O absorption, which respectively strengthens 
and weakens from luminosity class V to I (e.g. Aaronson 
et al. 1978; Kleinmann & Hall 1986; Origlia et al. 1993). 
The strength of various metal lines, e.g. the Na doublet 
at 2.20 /im and the Ca triplet at 2.26 /zm, also depend on 
the luminosity class. Ramirez et al. (1997) found that the 
quantity log[EW(CO)/(EW(Na) + EW(Ca))] is an excel- 
lent luminosity indicator, at least for K and M giants and 
dwarfs. Bessell et al. (1989a) identified the CN band at 
1.1 /im as a surface gravity sensitive feature. 

In this section, we summarise the results of a sys- 
tematic search for indices that separate the giants from 
supergiants of our sample. Our sample contains 10 su- 
pergiant observations: 3 non-variable supergiants, plus 3 
Semi-Regular variables of the SRc type, plus 3 Irregular 
variables of the Lc type, one of them observed twice. All 
these supergiants are of spectral type MO or later. The 
sample stars being cooler than those of Ramirez ct al. 
(1997) and our spectral resolution lower, the Cal and Nal 
features around 2.2 /im are seriously blended with molec- 
ular absorption features, leading us to redefine Na and 
Ca indices. The passbands and indices used are defined 
in Table ||. The passbands for CO measurements around 
2.3 /im are taken from Kleinmann & Hall (1986). F180W 
and F171W correspond to the Hubble Space Telescope 
narrow band NICMOS filters. They were initially designed 
to measure the strength of the 1.7 /im C2 absorption band 
in carbon stars (Thompson, private communication), but 
actually provide a good H 2 absorption index (see Langon 
et al. 1999 for applications to stellar population studies). 
Other common H2O indices (e.g. based on Wing filters at 
2.0 and 2.2 /im) would lead to qualitatively similar results. 
The selected filters are displayed on Fig. |l0| . 

The CO, H2O and metal line measurements are plotted 
in Fig(s). [ll] and |l2|. Supergiants systematically provide 
the highest metal line and CO indices of the sample of 
O-rich stars (carbon stars are discussed in Sect. 8). As 
expected, their mean CO and metal absorption is high and 
their mean H2O absorption low, compared to the mean 
properties of the M type variable giants of luminosity class 



Table 1. Filters used to define indices liable to separate 
giants and Miras from supergiants 



Filter 


Center 


Width 


Comments 


name 


(A) 


(A) 




F171W 


17150 


700 


HST NICMOS band 


F180W 


18050 


700 


HST NICMOS band, H 2 + C 


xNaci 


21920 


20 


'continuum' # 1 for Na 


xNac2 


21960 


20 


'continuum' # 2 for Na 


xNa 


22036 


152 


Na I doublet 


xNac3 


22120 


20 


'continuum' # 3 for Na 


xNac4 


22168 


25 


'continuum' # 4 for Na 


xCaci 


22548 


75 


'continuum' # 1 for Ca 


xCa 


22630 


90 


Ca I triplet 


xCac2 


22703 


55 


'continuum' # 2 for Ca 


xCOc 


22899 


52 


'continuum' for CO 


xCOi 


22957 


52 


12 CO (2,0) band- head 


x0O 2 


23245 


54 


12 CO(3,l) band-head 



CO = -2.5 log[(xCOi+xC0 2 )/xCOc] -CO(Vega) 
Na = -2.5 log[xNa/(xNaci + . . . +xNac 4 )] -Na(Vega) 
Ca = —2.5 log[xCa/(xCaci+xCac2)] — Ca(Vega) 
H 2 = -2.5 log[F180W/F171W] -H 2 0(Vega) 




10000 15000 20000 

wavelength (A) 

Fig. 10. Filters used in Sect. 7. Same spectra as Fig. || 

III. The combined constraints CO > 0.51, H 2 < 0.15 
and (Na+Ca) > 0.09 isolate the supergiants with only 3 
intruders: one of the two observations of KVCar (SRb 
variable), one of the three observations of WWSco (Mira 
variable) and the observation of V774 Sgr (Lb variable) . 

The luminosity class sensitivity of the CO/ (Na+Ca) 
combination of Ramirez et al. (1997) indeed extends to 
giant and supergiant luminosities. However, since (at a 
resolution of 1100) the centers of the Na and Ca line blends 
lie typically 10 % below the so-called continuum and their 
range of variation is small compared to the range of vari- 
ation of CO absorption, the combined index carries infor- 
mation similar to the CO index itself. 
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Fig. 11. The H2O index versus the CO index (see Sect. 7). 
Same symbols as Fig. [I] except for M Semi-Regulars 
(empty circles: SRa; filled circles: SRb). Supergiants (tri- 
angles) are discriminated in this diagram 



The isolation of the supergiants appears particularly 
clearly when a CO measurement is combined with the 
gravity indicator of Bessell et al. (1989a), the CN index 
110-108 (Fig. [[3]). In the same way, TiO indices measuring 
the <t> bands at 1.1 or 1.2 fim are absent in supergiants, and 
provide an alternative segregation tool (Fig. |lj|). These 
TiO indices are useful to identify very cool AGB stars. 
The 108, 110 and TiO cj> filters are shown in Fig. p|. 
Figures [Ti- 



ll 



also indicate a tendency for SRa type 
LPVs to have a lower CO index than SRb type variables. 
According to Kerschbaum & Hron (1992, 1994), SRa type 
variables appear as intermediate objects between Miras 
and SRb variables in several aspects. This assertion seems 
supported by Fig(s). pd|-^4] concerning the CO absorp- 
tion band. We note nevertheless that only 2 SRa variables 
(adding up to 8 observations) and 4 SRb variables (10 ob- 
servations) have been observed. More spectroscopic data 
should be studied to confirm this effect. 



8. Chemical types 

The spectra of carbon stars have little in common with 
those of oxygen-rich giants appart from CO absorption 
and a red energy distribution (Fig. [l0|). While, accord- 
ing to both observations and model atmospheres, it might 
be possible to identify selected regions of the near-IR 
spectra of O-rich LPVs as quasi-continuum windows (e.g. 
1.04 /im, 2.15 fan), the whole optical and near-IR range 
of C-rich stars is shaped by molecular bands of carbon 
based molecules. In the wavelength range under consid- 
eration the C-rich spectra display much less temperature 
and pulsation phase dependence than O-rich spectra. 

Our sample includes 26 NIR spectra of 6 carbon-rich 
stars. Four are Semi- Regular LPVs (7 spectra of RUPup, 
6 of YHya, 3 of TCae, 1 of SCen). The two Mira-type 



Fig. 12. The (Na+Ca) index versus the CO index (see 
Sect. 7). Same symbols as Fig. [Tl] 
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Fig. 13. The CN absorption index 110-108 versus the CO 
index. Same symbols as Fig. [TTJ 
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CO 

Fig. 14. The absorption index of TiO 4> (Av=0) versus 
the CO index. Same symbols as Fig. |ll| 
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variables are the dust-enshrouded RLep and the SC star 
BH Cru. Although the sample may not be representa- 
tive of all C-rich stars, it allows us to demonstrate which 
spectrophotometric indices are useful in the separation of 
chemical subtypes. 

It is well known that stars can be to a certain point 
separated into carbon stars and oxygen stars using broad 
band colour-colour diagrams (e.g. Frogel et al. 1980, 1990; 
Feast 1994; Zijlstra et al. 1996; Loup et al. 1998). C-rich 
stars appear in general redder than O-rich stars in such 
diagrams, and as a consequence, an often heard statement 
is that carbon stars are the reddest giants. 

The present data shows that this is not in general 
true. It is based on Magellanic Cloud data and is likely 
to hold for metal-deficient environments, where few if any 
very cool O-rich AGB stars are expected to exist because 
dredge-up of carbon produces C-stars more efficiently, but 
it must be reconsidered in the solar neighbourhood or 
in more general surveys. Excluding objects of our sam- 
ple for which the presence of a circumstellar envelope is 
known, the O-rich LPV sequence reaches much redder val- 
ues than the C-rich sequence in most colours (e.g. V— I, 
V— K, R— K, I— K). The colour distribution is similar for 
both chemical types in H— K. 

Most pure optical/NIR colours are thus insufficient to 
enable an unambiguous discrimination if used alone, and 
colours strongly affected by the specific molecular features 
of carbon stars are needed. Zijlstra et al. (1996) suggest 
the infrared colour K— [12], in which C-stars appear rel- 
atively red due to HCN and possibly C2H2 absorption 
in the 12-14/im range (Hron et al. 1998). J— K also is 
a useful indicator. Due to CN and C2 absorption in the 
J band, the C-stars appear systematically redder in this 
colour than O-rich stars with, for instance, identical I— K 
values (Fig. 15). J— K > 1.5 or 1.6 has been used for the 



separation of C-stars in LMC surveys (Loup et al. 1998; 
Wood et al. 1999). O-rich spectra with higher J— K values 
exist among the stars with no obvious dust shells of our 
sample, but they are rare. They generally correspond to 
Mira-type variables observed at minimum light. The J— K 
criterion remains applicable at solar metallicities, but it 
misses the warmer of the C-stars and only holds if M- 
stars with strong mass loss (and circumstellar emission) 
have previously been excluded. The narrow-band molecu- 
lar indices discussed below avoid confusion. 

The most obvious features identifying C-star spectra 
are those of CN and C2 bands. Indices for CN measure- 
ments include the Wing pairs 110—108 or 110—123, as well 
as the more selective narrow band index based on the first 
2 filters in Table |[ Indices sensitive to the steep Ballik- 
Ramsay C2 bandhead at 1.77 fim (Ballik & Ramsay 1963; 
Goorvitch 1990) include the HST NICMOS intermediate 
band index F180W— F171W and the narrow band index 
based on the nearly adjacent last 2 filters in Table |[ The 
filters are shown in Fig. [L6l 



* C Miras 



C Semi-Regulc 



++++ 



Fig. 15. J— K versus I— K. Carbon-rich Miras are rep- 
resented by asterisks and carbon-rich Semi-Regulars by 
filled circles. Otherwise, same symbols as Fig. |^ 

Table 2. Filters used to define indices which measure car- 
bon species 



Filter 


Center 


Width 


Comments 


name 


(A) 


(A) 




xCNcont 


10840 


30 


'continuum' for CN 


xCN 


10885 


40 


CN A-X Av=0 transition 


xC2,cont 


17570 


100 


'continuum' for C2 


xC 2 


17750 


140 


Ballik-Ramsay band 



In Fig. 0, the NICMOS index F180W-F171W is plot- 
ted against the narrow band C2 index. The good corre- 
lation among each sequence comes from the overlap be- 
tween the passbands used for both indices. However, only 
in the narrow band index are the two filters close enough 
not to be sensitive to the broad absorption band of H2O. 
While the NICMOS colour is able to identify LPVs in- 
dependently of their chemical type, it is not capable of 
separating C-rich objects from O-rich ones (Langon et al. 
1999). 

Figure [l8| shows the CN index against the C2 index. 
Once again, the oxygen-rich stars are clearly separated 
from carbon-rich stars as expected. The O-rich stars are 
now grouped together instead of forming a sequence. 

Another noticeable feature in Fig(s). [IT] and |l8| is the 
apparent segregation between carbon Miras and carbon 
Semi-Regular variables. If, on average, Mira-type pulsa- 
tion is reached at a later evolutionary stage on the AGB 
than Semi-Regular pulsation, as suggested by recent stud- 
ies of period-luminosity sequences in the LMC (Wood et 
al. 1999), and if the C2 band strength depends primarily 
on the C/O ratio, which on average increases with time 
along the AGB, then one expects to find the C-rich Miras 
at higher C2 indices than the Semi- Regulars. This sim- 
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Fig. 16. Filters used in Sect. 8. Same spectra as Fig. || 
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Fig. 17. The NICMOS index F180W-F171W versus the 
C2 index. Same symbols as Fig. and Fig. n9 
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Fig. 18. The CN index versus the C2 index. Same symbols 
as Fig. and Fig. R~5| 
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Fig. 19. The 110-123 index versus the 236-220 index. 
Same symbols as Fig. [| and Fig. [If]. The six observations 
of the C-Mira RLep and the three observations of the SC- 
Mira BH Cru are indicated. The arrow gives the extinction 
vector for an Eb-v equal to 1 



pic global scenario cannot be applied to samples as small 
as this one. In particular, the SC star BHCru pulsates 
like a Mira variable although its C/O ratio is close to 1. 
Larger and more representative samples are required to 
study systematic effects of the pulsation type. 

The CO indices measured for carbon stars cover the 
full range observed for O-rich objects, including supergiant 
stars. In particular, the 2.3 /zm CO absorption of the SC 
star BH Cru compares to the strongest band observed in 
our supergiant sample. This must be taken as a warning 
against the exclusive use of the CO index in searches for 
red supergiant stars (although the rarity of SC type stars 
prevents this finding from being of major concern in the 
study of integrated spectra of stellar populations). 

As already mentioned in the case of C2, the compar- 
ison of the CO index defined in Sect. 7 with the larger 
baseline index 236—220 illustrates the advantages of using 
narrow filters located very close to the feature of interest 
(Fig. |l^): the extremely low 236—220 index of RLep is 
merely due to its extremely red energy distribution; the 
systematic shift in CO between the 236—220 values of C- 
rich and O-rich stars is explained by additional molecular 
absorption in the 220 band; the shift in both CO and CN 
between Bulge and solar neighbourhood objects is mainly 
the effect of extinction. Colour corrections are required to 
circumvent ambiguities in the large baseline indices (Fro- 
gel & Whitford 1987). 

9. 12 C/ 13 C ratio 

The 12 C/ 13 C ratio in red giants is an important indicator 
of nuclear processing that occurs inside a star. The first 
dredge-up occurs when the star ascends the first giant 
branch: the number ratio 12 C/ 13 C drops from its initial 
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Fig. 20. Spectra exhibiting 12 CO and 13 CO lines with de- 
creasing strength (solid line: BH Cru, an SC Mira; dashed 
line: RS Hya, an M Mira; dotted line: S Car, an M Mira) 



value of ~ 90 to between 18 and 26 (Charbonnel 1994). 
The intermediate- mass stars (> 3-4 M.q) undergo a sec- 
ond dredge-up at the beginning of the AGB: the 12 C/ 13 C 
ratio decreases a little bit more for these stars. Then, the 
12 C/ 13 C ratio increases along the Thermally-Pulsating 
Asymptotic Giant Branch as the third dredge-up enriches 
the stellar surface with freshly synthesized elements (e.g. 
Forestini & Charbonnel 1997). However, when Hot Bot- 
tom Burning occurs, this ratio drops again to small values 
(Boothroyd et al. 1993). The 12 C/ 13 C isotope ratio is thus 
likely to exhibit a large spread of values for a given sample 
of cool giants. Our set of spectra enables us to check for 
an evolution of the 12 C/ 13 C ratio. 

Figure ^ shows three examples of star which exhibit 
12 CO and 13 CO lines with decreasing strength (solid line: 
SC4.5-SC7 Mira BH Cru; dashed line: M6 Mira RSHya; 
dotted line: K5-M6 Mira SCar). Positions of the CO 
bandheads are shown. The spectra are normalized to a 
common flux (F A ) around 22 900 A. 

A semi-quantitative approach using photometric in- 
" A 12 CO / 13 CO index is corn- 



dices is illustrated in Fig. 21. 

puted from the fluxes in rectangular 80 A wide filters, re- 
spectively centered in the four 12 CO and the three 13 CO 
bandheads shown in Fig. It is expected to vary with 
the 12 C/ 13 C ratio. It is plotted in Fig. pl| against the C 2 
index (Table ^) which separates C-rich from O-rich stars. 
Interstellar extinction has a negligible effect on the indices 
of Fig. 0. 

The C-rich stars exhibit the larger 12 CO / 13 CO indices: 
as expected, the 12 C/ 13 C ratio increases as the star un- 
dergoes third dredge-up that turns it into a carbon star. 
Supergiants (triangles) have also large 12 CO/ 13 CO indices 
compared to giants and LPVs. It may not be due to in- 
trinsic large values of 12 C/ 13 C. Indeed, predicted 12 C/ 13 C 
ratios from computations are weakly dependent on stellar 
mass at M. > 2 Mq (El Eid 1994): supergiants are not 



Fig. 21. The 12 CO/ 13 CO index as a function of the C 2 
index. Same symbols as Fig. [| and Fig. p"5| 

expected to exhibit large 12 C/ 13 C ratios. On the contrary, 
supergiants have in most cases low observed 12 C/ 13 C 
ratios (a Ori and a Sco for instance, Harris & Lambert 
1984). More probably, the large 12 CO/ 13 CO are due to 
radiative transfer effects on the CO bands (luminosity, at- 
mospheric extension, line blends with other species). 

Oxygen-rich LPVs exhibit a large spread of 
12 CO/ 13 CO indices, sometimes with values compa- 
rable to C-rich stars. Model atmospheres would be 
necessary to carefully determine the 12 C/ 13 C ratios 
of our sample stars and explain some of the trends in 
Fig. 21, and compare to predicted computations from 



evolutionary models. 



10. Conclusion 

We have explored the large sample of optical and near- 
infrared spectra of Langon & Wood (1998 and in prepa- 
ration), and derived useful tools to connect spectropho- 
tometric signatures with fundamental properties of cool 
oxygen- and carbon-rich LPVs and supergiants. 

Temperatures of oxygen-rich stars are assigned on 
the basis of synthetic spectra, computed with up-to-date 
oxygen-rich model atmosphere grids (Plez 1992; Plez, pri- 
vate communication) and adjusted to the optical data. 
They range from 2500 K to 4000 K. It is recalled that 
the relationship between this semi-empirical temperature 
and effective temperature is expected to depend on pulsa- 
tion properties as well as on the more usual physical at- 
tributes of a star (metallicity, surface gravity) . While sev- 
eral optical colours appear to be rather good indicators 
of the assigned temperature (e.g. 710—888), only a few 
purely near-infrared temperature indicators were found. 
The Wing colour 105—104 and the broad-band colours I— J 
and I— K can be used as a first approximation. The lack 
of better near-infrared temperature indicators is unfor- 
tunate for extragalactic survey purposes. V— K, although 
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not a purely NIR colour (and thus also more sensitive to 
extinction), proves to be a very convenient indicator of 
the optically assigned temperature over the whole range 
exhibited by the sample stars. 

The large spectral region covered by the observations 
enables us to derive a relation between the bolometric cor- 
rection BCi and the broad-band colour I— J for oxygen- 
rich cool stars. The standard deviation is only of 0.09 mag. 
It can be applied as a first approximation to red giants, su- 
pergiants and LPVs. This relation is of particular interest 
for ongoing infrared sky surveys. 

Interstellar extinction towards most stars of the so- 
lar neighbourhood sample appears to be small. Only four 
local stars (IOVir, WXPsc, RLep and IRC -20427) are 
clearly dust enshrouded. Even the latter display deep CO 
bands at 2.3 /xm, arguing against a significant contamina- 
tion by circumstellar dust emission in the near-IR. 

Deep H 2 features present in the near-IR spectra are 
a strong indication of pulsation. Nevertheless, the absence 
of water vapour does not exclude that the observed star 
is an LPV. 

Luminosity class I stars can be efficiently distinguished 
from luminosity class III stars with appropriate combina- 
tions of indices based on CO, H 2 0, Na, Ca, CN or TiO 
lines. SRa type LPVs show a tendency for lower CO in- 
dices than SRb variables. 

In order to separate carbon-rich LPVs from oxygen- 
rich stars, colours strongly affected by the specific molec- 
ular features of carbon stars are needed. Selected broad- 
band colours (e.g. J— K) are useful indicators, but narrow- 
band molecular indices based on CN and C 2 bands are by 
far more efficient. 

A 12 CO/ 13 CO index has been computed from the em- 
pirical spectra. As expected if the 12 CO / 13 CO index actu- 
ally varies with the 12 C/ 13 C ratio, the C-rich stars exhibit 
on average larger index values than the O-rich stars. The 
supergiants of the sample also exhibit large 12 CO/ 13 CO 
indices. This behaviour is opposite to what is expected 
from evolution models. This might be due to radiative 
transfer effects on the CO bands (luminosity, atmospheric 
extension, line blends with other species). 

The different indices and relations discussed are useful 
in the interpretation of integrated NIR spectra of stel- 
lar populations and of the considerable amount of data 
provided by ongoing infrared sky surveys, as well as in 
infering stellar properties of individual cool giants and su- 
pergiants. However, in order to resolve remaining ambigui- 
ties and to gain quantitative information, the dependence 
of the observable indices on stellar parameters, surface 
abundances, environement, etc, must be more extensively 
investigated. Only comparisons to models will eventually 
permit to complete this work. Extreme difficulties result 
from the fact that many complex aspects of LPV mod- 
elling (oxygen- rich/carbon-rich transition, pulsational dis- 
tinction between Mira and Semi- Regular variables, effects 
of pulsation on the atmosphere, non-grey emergent spec- 



tra, etc.) should be taken into account simultaneously. 
Some recent and important progress allows us to hope 
that this formidable task might be carried out in a future 
not completely out of reach. 
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